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Abstract 

Various approaches have been developed over the last four decades to characterize 
the magnitude of nonlinear interactions between triads of internal waves in stratified 
oceanic flows. The present manuscript compares some of these approaches and their 
predictions for internal wave nonlinearity parameter and Boltzman rate. We demon- 
strate that, for resonant triads in the limit of long internal waves in hydrostatic balance 
and in the absence of rotation, these various approaches predict equivalent rates of en- 
ergy transfer between waves. However, with the inclusion of background rotation and 
off-resonant interactions, these approaches lead to qualitatively different predictions. 
In particular, a noncanonical approach in Lagrangian coordinates leads to higher levels 
of nonlinearity at high frequencies and large wavenumbers than a canonical approach 
in isopycnal coordinates. 

1 Introduction 

Wave-wave interactions in stratified oceanic flows have been a subject of intensive research 
in the last four decades. Of particular importance is the existence of a "universal" internal- 
wave spectrum, the Garrett and Munk spectrum. It is generally perceived that the existence 
of a universal spectrum is, at least in part and perhaps even primarily, the result of nonlinear 
interactions of waves with different wavenumbers. Due to the quadratic nonlinearity of the 
underlying primitive equations and the fact that the linear internal-wave dispersion relation 
can satisfy a three-wave resonance condition, waves interact in triads. Therefore the question 
arises: how strongly do waves with wavenumber p interact with wavevectors pi and P2, where 
p = pi + p 2 ? What are the oceanographic consequences of this interaction? 
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Vari ous approaches have been develop ed 



actions 


Hasselmann . 


1966; 


Kenvon . 


1966 




L968: 


Vic Comas. 


1975; 


Miiller anc 


Olbers. 


1975 




Olbers, 


1974. 


1976 


; Pelinovskv and Raevskv 


. 1977 


McComas and Bretherton. 


1977 


Pomphrev et al. 


1980; 


Voronovich. 


1979; 


Milder. 


19821: Caillol and Zeitlin. 2000: 


Lvov and Tabakl. 


2001. 


2004) 



(see Table [U for a summary of the major distinctions.) 

All these approaches represent various attempts to derive a closed equation representing 
the slow time evol ution of the wave field' s wave action spectrum. Such an equation is called 



a kinetic equation ( iZakharov et all Il992l ). 



Table 1: A list of v arious kinetic equation s. 



Results from 
(|1980D are reviewed 



McComas and Bretherton! (119771); IPomphrev et al. 

( 1986 ). who state that Olbersl ( 1976 ). McComas and Bretherton ( 1977) and an unspec ified 



Olbers ( 


1976 


); 


in Miiller et al. 



Eulerian representation are consistent on the resonant ma nifold. Pomphrev et al. ( 1980h uti- 



lizes Langevin techniques to assess nonlinear transports. Miiller et al\ (119861 ) characterizes 
those Langevin results as being mutu a lly consistent with the direct eva lu ations of kinetic 



equations presented in lOlberd (119761); iMcComas and Bretherton (119771 ). iKenyonl (119681 ) 
states (without detail) that IKenyonl (119661 ) and iHasselmannl (119661 ) give numerically similar 
results. A formulation in terms of discrete modes will typically permit an arbitrary buoyancy 
profil e, but obtaining results requires specification of the profile. Of the discrete formula- 
tions, IPomphrev et al. (jl980h use an exponential profile and the others assume a constant 
stratification rate. The kinetic equations marked by ' are investigated in £JH while kinetic 
equations marked by * are investigated further in £j5l 



source 


coordinate 


vertical 


rotation 


hydro- 


special 




system 


structure 




static 




Hasselmann (1966) 


Lagrangian 


discrete 


no 


no 




Kenvon (1966. 1968) 


Eulerian 


discrete 


no 


no 


non-Hamiltonian 


Miiller and Olbers (1975) tJ 


Lagrangian 


cont. 


yes 


no 




McComas (1975. 1977) 


Lagrangian 


cont. 


yes 


yes 




Pelinovskv and Raevskv (1977) 


Lagrangian 


cont. 


no 


no 


Clebsh 


Voronovich (1979)^ 


Eulerian 


cont. 


no 


no 


Clebsh 


Pomphrev et al. (1980) 


Lagrangian 


discrete 


yes 


no 


Langevin 


Milder (1982) 


Isopycnal 


n/a 


no 


no 




Caillol and Zeitlin (2000) f 


Eulerian 


cont. 


no 


no 


non-Hamiltonian 


Lvov and Tabak (2001) t 


Isopycnal 


cont. 


no 


yes 


canonical 


Lvov and Tabak (2004) 1 


Isopycnal 


cont. 


yes 


yes 


canonical 



In this manuscript we concentrate on four of these different versions of the internal-wave 
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kinetic equation: 



a non canonical description using Lagrangian coordinates (lOlbera . 1 1 9 74l . 1 1 9761 : iMiiller and Olbers 



1975|), 



a canonical Hamil tonian description using Clebsh variables in Eulerian coordinates 
(jVoronovichl . Il979l ). 



a dynamical derivatio n of a kinetic equat i on wi thout use of Hamiltonian formalisms in 
Eulerian coordinates (jCaillol and Zeitlinl . |2000| ) , 



a can onical Hamiltonian description in isopycnal coordinates (jLvov and Tabakl . l200ll . 



20041 ). 



Our intent is to compare these approaches, and in particular, compare the predictions for 
the wavenumber- dependent characteristic nonlinear time scale of the Garrett and Munk wave 
action spectrum. To achieve this goal, we give necessary background in Section [2], briefly 
review approaches of Table [1] in Section ([3]), and then we demonstrate in Section H] that, 
under assumption of hydrostatic balance and under the assumption of resonant wave-wave 
interactions, the interaction matrices associated with the listed approaches are equivalent. 
While one, with sufficient experience, might regard this as an intuitive statement, it is far 
from trivial. We will then demonstrate in Section [5] that, if the assumption of resonant wave- 
wave interactions is relaxed, both quantitatively and qualitatively different transfer rates are 
predicted. In particular, we show that the Boltzman rate e p , defined below in (j2J) is, in fact, 
representation dependent for near-resonant interactions. 

We have not, at this time, achieved a detailed mathematical understanding of how these 
differences arise and consequently do not digress into a detailed discussion of why, for exam- 
ple, the radius of convergence of two consecutive series expansions in one coordinate system 
differs so dramatically from a single series expansion in a different coordinate system. We 
conclude in Section [6j 



2 Background 

A kinetic equation is a closed equation for the time evolution of the wave action spectrum 
in a system of weakly interacting waves. It is usually derived as a central result of wave tur- 
bulence theory. The concepts of wave turbulence theory provide a fairly general framework 
for studying the statistical steady states in a large class of weakly interacting and weakly 
nonl inear many-body or m any-wave systems. In its essence, classical weak turbulence the- 



ory (jZakharov et all Il992l ) is a perturbation expansion in the amplitude of the nonlinearity, 
yielding, at the leading order, linear waves, with amplitudes slowly modulated at higher or- 
ders by resonant nonlinear interactions. This modulation leads to a resonant redistribution 
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of the spectral energy density among space- and time-scales, and is described by a kinetic 
equation. 

Typical assumptions needed for the derivation of kinetic equations are: 

• Weak nonlinearity, 

• Gaussian statistics of the interacting wave field in wavenumber space and 

• Resonant wave-wave interactions. 

The derivation of the kinetic equation for general nonlinear systems is well studied and 
under stood, and thus will not be repeated here. Three wave kinetic eq uations take the 
form (jZakharov et all Il992l ; iLvov and Nazarenko 1 . 120041 : iLvov et all Il997l ) : 



r 

£ = 4tt / iV^pJ 2 f pl2 5 P - Pl - P2 5(u p - u Pl - Up 2 )dp 12 
-4vr J \V£l p \ 2 f 12p <5 Pl _p 2 _p50 Pl - u P2 - u p ) dp 12 



-Air 



I ypz 1 2 



/2pl 8 f 



p,pi I J ^P 1 "P2—P—P1 

with /. 



5(u. 



pi 



10. 



Pi 



dp 



12 , 



pl2 ^lp\^p2 Tlpijlpi Tlp2*) ■ 



Here n p = n(p) is a three-dimensional wave action spectrum (spectral energy density divided 
by frequency) and the interacting wavevectors p, pi and p 2 are given by 



V 



[k,m), 



i.e. k is the horizontal part of p and m is its vertical component. We assume the wavevectors 
are signed variables and wave frequencies u p are restricted to be positive. The magnitude 
of wave-wave interactions V p p pi is a matrix representation of the coupling between triad 
members. It serves as a multiplier in the nonlinear convolution term in what is now commonly 
called the Zakharov equation - equation in the Fourier space for the waves field variable. 
This is also an expression that multiplies the cubic convolution term in the three-wave 
Hamiltonian. 

For internal waves in the ocean such kinetic equation was derived by the approaches in 
Table [TJ The development of a kinetic equation is facilitated by transforming to canonical 
coordinates in a Hamiltonian framework, for which one can demonstrate that the symme- 
tries and hence conservation principles of the original equ ation set in the spatial /temporal 
domains have been preserved in the spectral domain (e.g. Zakharov et all Il992 ). Finding 
canonical coordinates, however, can be highly nontrivial. Transformations to canonical coor - 
dinates have been found usin g Clebsch variables in Euleria n coordinates (jVoronovichl . Il979l ) 
and in isopycnal coordinates (jLvov and Tabaki . l200ll . 120041 ) . Kinetic equations in Lagrangian 
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coordinates start by averaging the Lagrangian of the stratified fluid (averaging of the vari- 
ational principle) and then transform the Lagrangian to the Hamiltonian. The Lagrangian 
coordinate kinetic equations considered here are noncanonical. We also note that it is pos- 
sible to obtain a kinetic equation directly from the dynamical equations of motion, without 
the us e of the Hamiltonian structure. Such an approach was executed by ICaillol and Zeitlin 
(120001 ). The conservation properties of non-canonical and non-Hamiltonian representations 
are not guaranteed unless explicitly demonstrated. The issue of conservation properties is 
greatly compounded for non-canonical and non-Hamiltonian representations off the resonant 
manifold. 

A typical restriction is to exclude interactions with potential vo rticity carrying m e mbers 
of t he fluid dynamical system , for which we refer the reader to iLelong and Riley! (1l99ll ) 
and ICaillol and Zeitlin! (120001 ) . Even in these extended analyses a plane wave formulation is 
assumed t hat eliminates the pot e ntial vorticity assoc iated with a slowly varying wave-packet 
structure dBiihler and Mclntvrel . 120051 : IPolzinl . l2008h . 

Note that the kinetic equation allows us to numerically estimate the life time of any 
given spectrum. In particular, we can define a wavenumber dependent nonlinear time scale 
proportional to the inverse Boltzman rate: 



.NL 



hr 



(2) 



This time scale characterizes the net rate at which the spectrum changes and can be directly 
calculated from the kinetic equation. 

One can also define the characteristic linear time scale, 



2ir/u p . 



The non-dimensional ratio of these time scales can characterize the level of nonlinearity in 
the nonlinear system: 



-NL 



2-nhr 



n p uj p 



(3) 



We refer to (J3]) as the nonlinearity parameter. 

The nonlinear parameter serves as a low order consistency check for the various kinetic 
equation derivations. An 0(1) value of e p implies that the derivation of the kinetic equation is 
internally inconsistent. The Boltzman rate represents the net rate of transfer for wavenumber 
p and is an appropriate measure of nonlinearity for smooth, isotropic and homogeneous 
spectra. The i ndividual rates of transfer into and out of p maybe significantly larger for 
spectral spikes ( IMuller et all Il986l ) and potentially for smooth, homogeneous but anisotropic 



spectra. Estimates of the Boltzman rate and e p require integration of Eq. ([I]). In this 
manuscript such integration is performed numerically. 
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In this paper we concentrate on four a p proaches, namely IMiiller and Olbersl (Il975l ): 
Voronovichl (Il979h : ICaillol and Zeitlinl fl2000h : iLvov and Ta~bakf fl200ll . l20Q4h . We show that 



on the resonant manifold they produce equivalent results. 

Resonant interaction approximation is self-consistent for small level of nonlinearities. 
However, as the nonlinearity parameter increase, near-resonant interactions start to play a 
role. 

For realistic estimates the effects of rotation must be included, and this restricts our 
investigations to two approaches tha t allow inclusion of bac kgro und rotations. Therefore , 
we concentrate in more details on the lLvov and Tabaki (12004 ) and lMiiller and Olbersl (119751 ) 
representations. 

We show that for the near-resonant interactions, these two approaches returns qualita- 
tively different predictions for transfer rates. This is the main physical result of the present 
paper. 

There is a multitude of reasons for possible differences. First and foremost, we view the 
distinction between Lagrangian, isopycnal and Eulerian coordinates as the most dynamically 
significant difference. The use of a Lagrangian coordinate system requires an expansion in 
powers of small fluid parcel displacements in addition to an assumption of weak nonlinearity, 
whereas formulations in isopycnal or Eulerian coordinates require only an assumption of weak 
nonlinearity. An issue with extant Lagrangian coordinate representations is that the small 
amplitude assumption represents an unconstrained approximation whose domain of val idity 



vis-a-vis the weak interaction approximation is not well defined, (IMiiller et all [1986). A 



subsidiary issue is that the use of a Lagrangian coordinate system places the nonlinearity in 
the incompressibi lity constra i nt, an d a single plane wave is not an exact solution of the equa- 
tions of motion, (jSandersonl . 119851 ). Similarly, a single plane wave also does not constitute 
a solution to the isopycnal equations of motion. In Eulerian coordinates the nonlinearity is 
advective and a single plane wave is an exact solution of the equations of motion. On the 
other hand, it is a robust observational fact that Eulerian frequency spectra at high vertical 
wavenumber are contaminated by vertical Doppler shifting: near-inertial frequency energy 
is Doppler shifted to higher frequency at approximately the same ver tical wavelength. Use 
of an isopycnal coordinate system considerably reduces this artifact, (jSherman and Pinkell . 



19911 ) . Thus differences in the approaches may represent physical effects rather than techni- 



cal issues such as the prop er implementation of a potential vorticity conservation statement 
( ICaiUol and Zeitlinl . [iooof l . 

We emphasize that our intent is to estimate transport rates for various approaches within 
a common framework and to compare those results. Our goal is a qualitative physical expla- 
nation of the possible reasons for the similarities and differences rather than a quantitative 
analytical explanation of how those differences arise. 
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3 Various Approaches 



In this section we list the approaches that we use. We do so for completion and to transfer 
everything to a uniform notation. Our attention is restricted to the hydrostatic balance case, 
for which 

|fe| <C \m\ . (4) 

A minor detail is that the linear frequency has different algebraic representations in isopy- 
cnal and Cartesian coordinates. The Cartesian vertical wavenumber, k z , and the density 
wavenumber, m, are related as m = —g/(poN 2 )k z where g is gravity, p is density with refer- 
ence value po, N is the buoyancy (Brunt-Vaisala) frequency and / is the Coriolis frequency. 
In isopycnal coordinates the dispersion relation is given by, 



plN 2 m? 



(5) 



In Cartesian coordinates, 



«(p) = \ f 2 + N 2 



\k\ 2 
p 

h z 



(6) 



In the limit of / = these dispersion relations assume the form 

\k\ \k\ 

Up OC r OC — — - 

\m\ \k z \ 



(7) 



3.1 Kenyon and Hasselmann 

The first kinetic equatio n s for wave-w a ve int eract ions in a conti nu ously st r atified ocean 
appear in Kenyonl (Il966 ). Hasselmann ( 1966) and Kenvon ( 1968 ). Kenyon ( 1968 ) states 



(without detail) that Kenyonl ( 1966[) and Hasselmann ( 19661 ) give numerically similar results 



We have found that iKenyonl (119661 ) differs from the four approaches examined below on one 
of the resonant manifolds, but have not pursu ed the questio n further. It is possible this 
difference results from a typographical error in IKenyonl (119661 ). We have not rederived this 
non-Hamiltonian representation and thus exclude it from this study. 



3.2 Miiller and Olbers 



Matrix elements derived in lOlbersI (Il974t ) are give n by W^ T . MO \ 2 = T + J^ and l^ p M °l 2 = 
T~/ (47r). We extracted T ± from the Appendix of Miiller and Olbers ( 19751 ). In our notation, 
in the hydrostatic balance approximation, their matrix elements are given by 
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\v p 

1 Pi ,P2 



MO|2 



(iVo 2 - P) 
32p 



■UJLO\<jJ<i 





k\ 


fel 


&2 




UJUJiUJ2\p\ 


Pll 


|P2| 



-mi 



ki-k'2—ifk2-k^/uji 



+ m 2 -m 2 



ki-k'2— ifhi-hi- Iu2 . 

fc| hm l 



m 



k2-k+ifk2-k ± /uj2 1 A f k2-k—ifh-hn/u> , 

-m 2 J k 2 — — + m) -m J k2 2 + m 2 



mi 



-m 



kki —ifk-k^ /lu 



+ m i ) —mi 



+ m 



m 2 



Taking a / = limit we get: 



I pi ,P2 



MOi2 



oc 



k\ 


Ifcll 




1 


\ m \ 


\mm\m2\ 





1 ^ m 2 fc ■ ^2 

mi \ 



Ifc9 12 



mfc 2 • fc 



+ m 2 M 

m 2 



m 2 fci • fc 2 
|fc 2 | 2 + 
mfei • A; 
Ifel 2 



mi 



mi 



m 1 k 2 ■ ki 

rriik ■ k\ 
Ifcil 2 " 



+ m 2 



m 



(9) 



3.3 Pelinovsky and Raevsky 



An important paper on internal waves is iPelinovsky and Raevskyl (119771 ). Clebsh variables 
are used to obtain the interaction matrix elements for both constant stratification rates, 
iV = const., and arbitrary buoyancy profiles, N = N(z) . Not much details are given, but 
there are some similarities in appearance with Voronovichl ( 1979 ). The most significant result 
is the identification of a scale invariant (non-rotating, hydrostatic) stationary state. It is 
stated in the p aper that th e ir ma trix elemen ts are equiv a lent to thos e derived in their citation 
[11], w hich is Brehovsky ( 1975 ). Because Brehovskyl ( 19751 ) and Pelinovsky and Raevsky 
( 119771 ) are in Russian and not generally available, we refrain from including them in this 
comparison. 



3.4 Voronovich 

Voronovich used Clebsh-like variables to derive the Hamiltonian for incompressible stratified 
flows in the ocean. It is probably the first canonical Hamiltonian structure derived for such 
kind of fl ows. A deta i led ex planation of Voronovich's method appears in section 7.1 of the 
textbook iMiropolskyl (119811 ) It is a straightforward task to write down the kinetic equation 
associated with this Hamiltonian structure. 
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We formulate the matrix elements for Voronovich's Hamiltonian using his formula (A.l). 
This formula is derived for general boundary conditions. To compare with other matrix 
elements of this paper, we assume a constant stratification profile and Fourier basis as the 
vertical structure function (f>(z). That allows us to solve for the matrix elements defined 
via Eq. (11) and above it in his paper. Then the convolutions of the basis functions give 
delta- functions in vertical wavenumbers. Vornovich's equation (A.l) transforms into: 



\V P 



V 1 2 I ^1 1 ^1 1 1 ^2 | 
OC -j r 

|mm 1 m 2 | 



-777. 



1 



\k to 



+777,1 

+m 2 



1 



fei mJ 



k 


k\ 777i 






k 


fei|m| 


\k\ 



k 2 \\m 2 \ 



k- 


k 2 


to 




k 





k ■ 

- + — 


k 2 \m 2 \ 




\k 2 \ 




k 2 \m 2 \ 




IN 


k 2 ■ k x 


777i|^ 


M J 



Ul\ + UJ 2 — UJ 

UJ 

Ul\ + uj 2 — UJ 

UJ 2 — UJ 



UJ 2 



(10) 



Note that Eq. (TTTJT) shares structural similarities with the interaction matrix elements in 
isopycnal coordinates, Eq. (1131) below. 



3.5 Milder 



An alternative Hamiltonian description was developed in iMilderl (119821 ) . in isopycnal coor- 
dinates without assuming a hydrostatic balance. The resulting Hamiltonian is an iterative 
expansion in powers of a small parameter, similar to the case of surface gravity waves. In 
principle, that approach may also be use d to ca l culate wave-wave interaction amplitudes. 
Since those calculations were not done in IMilderl (119821 ). we do not pursue the comparison 
further. 



3.6 Caillol and Zeitlin 



A non -Hamiltonian kinetic equation for internal waves was derived in ICaillol and Zeitlin 
(l2000h. Eq. (611. To make i t appear equivalent to more traditional form of kinetic equation, 



as m 
k -> 



Zakharov et al 



_ ( 119921 ). we m ake a change of var i ables I — > —I in the second line, and 
-k in the third line of (61) of ICaillol and ZeitlinI (l2000h . If we further assume that all 
spectra are sym metric, n{— p) = n(p), the n the kinetic equ a tion a s sumes traditiona l form 
as in Eq. 
2004T ). 



sec 



Miiller and Olbersl (119751 ): IZakharov et all (I1992T ): iLvov and Tabakl (12001 



The matrix elements according to ICaillol and ZeitlinI (120001 ) are shown as Xk t i tP and Y h 



in Eqs. (62) and (63), where \Vg 



CZ|2 



Pi-Pz 



Ap 1; p 2 ,p &nd 



yp 

P2,P 



CZ|2 



Y. 



P\,-P2,P' 



k,l,p 

In our notation 
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it reads 



\V P 

I pi ,P2 



cz | 2 oc (|fc|sgn(m) + |fc 1 |sgn(m 1 ) + |fc 2 |sgn(m 2 )) 2 j 

2 — |fci|sgn(m 1 )|A; 2 |sgn(m2) 



m 



mim 2 )' 



x 



fel 



m\ \rnx\ \m 2 \ |fe| |fci| \k 2 
fci 



-m 



2 ifc 2 r 2 



m^ — mim 2 



mi 



m 2 



(11) 



3.7 Isopycnal Hamiltonian 



Finally, in iLvov and Tabald (12004 ) the following wave- wave interaction matrix element was 
derived based on a canonical Hamiltonian formulation in isopycnal coordinates: 



IV? 



32^ 



k\k 2 



kk\ ■ k 2 \<jJ\<jJi k\k 2 ■ k fuJ^uJ k 2 k ■ k\ [uxJ[ 



OJ 



kok 



kk\ 



UJ 2 



f 2 k\k 2 ■ k — k\k ■ ki — k 2 k\ ■ k 2 



^/UJUJ\ijJ 2 



kk\k 2 



+ If 



kk\k 2 
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LOlL0 2 



■{k\ k 2 ) 



1 (e 2 - e) 



UJ 2 UJ 



UJiO\ 



ie-ki) 



(12) 



Lvov and Tabakl (120011) is a rotat ionless limi t of ILvov and Tabakl (12004f) . Taking the / — > 
limit, the lLvov and Tabakl (120041 ) reduces to ILvov and Tabakl (120011 ). and (IT2"|) reduces to 



\V P 

I Pl ,P2 



Hi2 



OC 



fe fei fe 2 





m 


fcjfci • k 2 J 


m\m 2 



+ lfellfc 2 ■ fel 



mi 



m 2 m 



+ |fe 2 |fe • fei' 



m 2 



mm i 



(13) 



Observe that in this form, these equations share structural similarities with Eq. ({TO]) . 

In this section we gave brief review of the various approaches that were developed for 
describing wave-wave interactions of internal waves in the ocean. While this review is neces- 
sarily brief, this is the first time all these papers are cited together by a single manuscript. 



4 Resonant wave-wave interactions 

How one can compare the function of two vectors p± and p 2 , and their sum or difference? First 
one realizes that out of 6 components of p\ and p 2 , only relative angles between wavevectors 
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enter into the equation for matrix elements. That is because the matrix elements depend on 
the inner products of wavevectors. The overall horizontal orientation of the wavevectors does 
not matter: relative angles can be determined from a triangle inequality and the magnitudes 
of the horizontal wavevectors k, k\ and &2- Thus the only needed components are |fc|, 
| fci|, Ifcal, m and ni\ (777.2 is computed from m and mi). Further note that in the / = and 
hydrostatic limit, all matrix elements become scale invariant functions. That is to say that it 
is sufficient to choose an arbitrary scalar value for |fe|, and m, since only |fci|/|fc|, |fc 2 |/|fc| and 
mi/m enter the expressions for matrix elements. We make the particular (arbitrary) choice 
that \k\ = m = 1 for the purpose of numerical evaluation, and thus the only independent 
variables to consider are |fei|, |&2| and mi. Finally, mi is determined from the resonance 
conditions, as explained in the next subsection below. As a result, we are left with a matrix 
element as a function of only two parameters, k\ and k^- This allows us to easily compare 
the values of matrix elements on the resonant manifold. 



4.1 Reduction to the Resonant Manifold 

When confined to the traditional form of the kinetic equation, wave-wave interactions (scat- 
tering) are constrained to the resonant manifolds defined by 

jp = Pl+P2 ,* fpi=P2+P v \P2=P + Pl 

a ) \ b ) \ C ) \ , ( 14 ) 

I LU = LO\ + LU2 I UJi = ix>2 + Ul \UJ2 = UJ + UJl 

To compare matrix elements on the resonant manifold we are going to use the above resonant 
conditions and the internal- wave dispersion relation (171). To determine vertical components 
mi and m2 of the interacting wavevectors, one has to solve the resulting quadratic equations. 
Without restricting generality we choose m > 0. There are two solutions for mi and m 2 
given below for each of the three resonance types described above. 
Resonances of type (TUb ) give 



mi 


m 1 
~ 2|fe| \ 


\k\ + 


I*i| 


+ \k 2 \ + 


m 2 


= 171 — 


mi. 






mi 


m 1 

~ 2\k\ \ 


f \k\ - 


Ifcil 


- Ifel - 


m 2 


= 171 — 


m\. 







V(|fc| + |fci| + |fc 2 |) 2 -4|fc||fci|) 

y/(\k\ ~ |fei| - |fc 2 |) 2 + 4|fc||fc 1 



(15a) 



(15b) 

Note that because of the symmetry, f)15al) translates to (|15bl) if wavenumbers 1 and 2 are 
exchanged. 

Resonances of type (fllb ) give 



2|fe| 1 |f " 

m + m 2 . 



(|fc| - |fei| - |fc 2 | + v/pi - |fci| - |fc 2 |) 2 + 4|fc||fc 2 | 



16a) 
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' m 2 = (\ k \ + \ k i\ " l fc 2| + y/(\k\ + |fei| - |fc 2 |) 2 + 4|fc||fc 2 

m 1 = m + m 2 . 

Resonances of type (lT4"b ) give 

^1 = " " N + v /(|fc|-|fc 1 |-|fc 2 |) 2 +4|fc||fc 1 

m 2 = m + To], 

mi = (|fe| - |fei| + |fca| + y/(\k\ - \k t \ + 1^1)2 + 4^1^! 
m 2 = m + mi. 



(16b) 



(17a) 



(17b) 



Because of the symmetries of the problem, (|16a[) is equivalent to (I17al) . and (116bj) is equivalent 
to (117bj) if wavenumbers 1 and 2 are exchanged. 



4.2 Comparison of matrix elements 

As explained above, we assume / = and hydrostatic balance. Such a choice makes the 
matrix elements to be scale-invariant functions that depend only upon |fei| and |fc 2 |. As a 
consequence of the triangle inequality we need to consider matrix elements only within a 
"kinematic box" defined by 

1 1 — |&2 1 1 < \k\ < | fci| + |fc 2 |. 

The matrix elements will have different values depending on the dimensions so that isopycnal 
and Eulerian approaches will give different values ([5])-([6]). To address this issue in the simplest 
possible way, we multiply each matrix element by a dimensional number chosen so that all 
matrix elements are equivalent for some specific wavevector. In particular, we choose the 
scaling constant so that |V(|fei| = 1, \kz\ = 1)| 2 = 1. This allows a transparent comparison 
without worrying about dimensional differences between various formulations. 



4.2.1 Resonances of the "sum" type ( 114a ) 

Figured] presents the values of the matrix element l^. P2 Ji5bl J 2 on ^ e resonan t sub-manifold 
given explicitly by (Il5bl) . All approaches give equivalent results. This is confirmed by 
plotting the relative ratio between these approaches, and it is given by numerical noise (not 
shown). The solution (I15al) gives the same matrix elements but with |fei| and |fc 2 | exchanged 
owing to their symmetries. 



4.2.2 Resonances of the "difference" type (114b) and ( 114c ) 

We then turn our attention to resonances of "difference" type (114b) for which (|14b ) could be 
obtained by symmetrical exchange of the indices. All the matrix elements | ^^JTKTl 1 2 on the 
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resonant sub-manifold (116aj) . are shown in Fig. [2j All the matrix elements are equivalent. The 
relative differences between different approaches are given by numerical noise (not shown). 
Finally, I^V|l6bjl 2 011 ^ ne resonan t sub-manifold (|16bl) are shown in Fig. [3) Again, all the 
matrix elements are equivalent. 

The solutions (I17al) and (|17bl) give the same matrix elements but with \k\\ and \k 
exchanged as the solutions ( 116ajl and ( 116bl) owing to their symmetries. 
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4.2.3 Special triads 



Three simple interaction mechanisms are identified by iMcComas and Brethertonl (119771 ) in 
the limit of an extreme scale separation. In this subsection we look in closer detail at 
these special limiting triads to confirm that all matrix elements are indeed asymptotically 
consistent. The limiting cases are: 

• the vertical backscattering of a high-frequency wave by a low frequency wave of twice 
the vertical wavenumber into a second high-frequency wave of oppositely signed vertical 
wavenumber. This type of scattering is called elastic scattering (ES). The solution ( 11 5 aft 
in the limit \ki\ ^ corresponds to this type of special triad. 

• The scattering of a high-frequency wave by a low-frequency, small-wavenumber wave 
into a second, nearly identical, high-frequency large- wavenumber wave. This type of 
scattering is called induced diffusion (ID). The solution (I15bj) in the limit that |fei| — ► 
corresponds to this type of special triad. 

• The decay of a low wavenumber wave into two high vertical wavenumber waves of ap- 
proximately one-half the frequency. This is called parametric subharmonic instability 
(PSI). The solution ( 116al) in the limit that |fei| — ► corresponds to this type of triad. 

To study the detailed behavior of the matrix elements in the special triad cases, we choose 
to present the matrix elements along a straight line defined by 

(|fc 1 |,|fe 2 |) = (e,e/3 + l)|fc|. 

This line originates from the corner of the kinematic box in Figs. [TH3]at (|fci|, |&2|) = (0, |fe|) 
and has a slope of 1/3. The slope of this line is arbitrary. We could have taken e/4 or 
e/2. The matrix elements here are shown as functions of e in Fig. HI We see that all four 
approaches are again equivalent on the resonant manifold for the case of special triads. 

In this section we demonstrated that all four approaches we considered produce equivalent 
results on the resonant manifold in the absence of background rotation. This statement is 
not trivial, given the different assumptions and coordinate systems that have been used for 
the various kinetic equation derivations. 
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5 Smearing of the resonance manifold and near-resonant 
Interactions 



5.1 Nonlinear frequency renormalization as a result of nonlinear 
wave- wave interactions 

Above we have compared the values of matrix element on the resonant manifold. The 
resonant interaction approximation is a mathematical simplification which reduces the com- 
plexity of the problem. In this subsection we examine transfers including near resonant 
interactions. Our interest in near-resonant interactions has significant physic al motivations . 



For example a major unresolved issue is the importance of Doppler shifting (IPolzin 1 . 12004 ). 



Of particular interest here is the variable effects of Doppler shifting in different coordinate 
systems. The resonant interaction approximation assumes, perforce, an expansion in terms 
of a non-advected wavefield, with dispersion relation given by Eq. (j5J) or Eq. (jBJ). In the 
limit of extreme time scale separation between high frequency waves and a low frequency 
background, one is tempted to replace the non-advected frequency by its Doppler shifted 
intrinsic frequency counterpart, u) — > u) — k ■ u, in which u and k are the frequency and 
wavevector of the high frequency wave and u is th e velocity field of t he low frequency 



wavefield. This is the genesis of the eikonal approach (jMuller et all , Il986l ) to internal wave- 
wave interactions. Then the resonant approximation is self-consistent for small values of 
nonlinearities. Indeed, change in the wave amplitude will be small, and the Doppler shift 
cancels from the frequency delta function. Yet, as nonlinearity increases, the near-resonant 
interactions become more and more pronounced, consequently the issue of Doppler shifting 
more and more important. Furthermore, near- resonant int eractions play a major role in 



numerical simulations on a discre te grids ( Lvov et al . 20061) , for time evolut ion of discrete 



ity waves (jJanssen 



2003 



systems (IGershgqrin et all 120070. in acoustic turbulence (jLvov et all Il997l). surface grav- 



Annenkov and Shrira. V.I 



20061 ) 



Yuen and Lakd . Il982l ) , and internal waves (I Voronovich et al.l . 12006 



To take into account the effects of near-resonant interactions self-consistently, the en- 
need to be "broadened" . The 



pi 



— LO. 



Vili 



ergy conserving delta-functions in Eq. (JTJ, 5(uj p — to. 
physical motivation for this broadening is the following: when the resonant kinetic equation 
is derived, it is assumed that the amplitude of each plane wave is constant in time, or, in 
other words, that the lifetime of single plane wave is infinite. Resulting kinetic equation, 
nevertheless, predicts that the amplitude of the wave do change. Consequently the wave life- 
time is finite. For small level of nonlinearity this distinction is not significant, and resonant 
kinetic equation constitutes a self-consistent description. For larger values of nonliterary 
this is no longer the case, and the wave lifetime is finite and amplitude changes need to be 
taken into account. Consequently interactions may not be strictly resonant. This statement 
also follows from the Fourier uncertainty principle. In other words, the waves with varying 



14 



amplitude can not be represented by a single Fourier component. This effect is larger for 
stronger level of nonlinearity parameter. 

Derivation of t he kinetic equation with a broadened delta function is given in details in 
( iLvov et all 119971 ). and is not going to be repeated here. The result is that 



dn p 



4 / \VI 12 



-4 
-4 



)dp 12 



pj.pjl /pl2 Op. 

\Vp°, Pl \ 2 hpi &P2-P-P1 £(, U P2 ~ u p ~ u pi) dpi2 , 



Here L is defined as 



£(Aw) 



kl2 



(19) 



where Tfc 12 is the total broadening of each particular resonance, and is given below. If the 
nonlinear frequency renormalization tends to zero, i.e. Tku — > 0, £ reduces to the delta 
function: 

lim C(Au) = tt5(Auj). 



Consequently, in the limit resonant in teractions (i.e. no broadening) (fl8l) reduces to Eq. 
. We have shown in ILvov et al\ (119971 ) that the broadening in Eq. (fl9l) is given by 



r fc i2 = 7p + 7pi + 7; 



1P2- 



(20) 



It means that the total resonance broadening is the sum of individual frequency broadening, 
and can be thus seen as the "triad interaction" frequency. 

The single frequency renormalization can be calculated self- consistently from 

7 P = 4 J \Vp° uP2 \ 2 (n pi + n P2 ) Vpi-p2 £(^p ~ w pi - u P2 )dp l2 

-4 J \Vp 2 ] P \ 2 (n P2 - n Pl )5 Pl - P2 - p C{u Pl - u P2 - u p ) dp 12 

-4 J \y p , 2 Pl \ 2 (n Pl - n P2 ) 5 P2 _ P _ P1 C(u P2 - u p - u pi ) dp 12 . 

(21) 

The interpretation of this formula is the following: nonlinear wave-wave interactions lead to 
the change of wave amplitude, which in turn makes the lifetime of the waves to be finite. 
This, in turn, makes the interactions to be near-resonant. 
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A self-consistent estimate of 7 P requires the iterative solution of (TIBl) and (|2Til over the en- 
tire field: the width of the resonance ff2Tl) depends on the lifetime of an individual wave [from 
018p ]. which in turn depends on the width of the resonance fl20|) . This numerically intensive 
computation is beyond the scope of this manuscript. Instead, we make the uncontrolled 
approximation that: 

7 P = 5u p . (22) 

We note that this choice is made for illustration purposes only, we certainly do not 
claim that it represents a self consistent choice. Below, we will take 5 to be 10~ 2 and 10~ 3 . 
These values are rather small, therefore we remain in the closest proximity to the resonant 
interactions. 



5.2 Characteristic nonlinear time scale of the Garrett and Munk 
Spectrum 

Estimates of near- resonant transfers are obtained by assuming horizontal isotropy integrating 
ffT8]) over horizontal azimuth: 

= / ^^'^' P2 ' 2 6p - pi ~ P2 C<yUJp ~ Upi ~ LJ P2) dk i2dm 1 
-An [ ^^{V^pl 2 f 12p 5 Pl - P2 - p C(uj pi - lu P2 - uj p ) dk 12 dmi 

J <->12p 

/k k 
P2 — uj p — ui Pl ) dkudmi , (23) 
>->2pl 

where S p ± 2 is the area of the triangle k = k\ + We numerically integrated ( 1231) for 
p's which have frequencies from f to N [specifically (33/32, 17/16, 9/8, 5/4, 3/2, 2, 4, 
...)/] and vertical wavenumbers from 47r/(26) to 2007r/(26) ([2,4,6, ... 98] n/b) . The limits 
of integration are restricted by horizontal wavenumbers from 27r/10 5 to 27r/10 meters -1 , 
vertical wavenumbers from 27r/(26) to 27r/10 meters -1 , and frequencies from / to N. The 
integrals over k\ and k 2 are obtained in the kinematic box in k\ — k 2 space. The grids in the 
k\ — k 2 domain have 2 17 points that are distributed heavily around the corner of the kinematic 
box. The integral over m\ is obtained with 2 13 grid points, which are also distributed heavily 
for the small vertical wavenumbers whose absolute values are less than 5m, where m is the 
vertical wavenumber. 

Below we calculate the nonlinear time scale ([2]) and nonlinearity parameter ([3]). To 
calculate this parameter, we need to choose a form of spectral energy density of internal 
waves. We we utilize the Garrett and Munk spectrum greed-upon representation of 

the internal waves: 
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Here the reference wavenumber is given by 



m* = 7rj'*/6, 



(25) 



in which the variable j represents the mode number of an ocean with an exponential buoyancy 
frequency profile having a scale height of b. 
We choose the following set of parameters: 

• b = 1300 m in the GM model 



The total energy is set as: 



E = 30 x 10 -4 m 2 s" 



• Inertial frequency is given by / = 10 4 rad/sec, and buoyancy frequency is given by 
iVo = 5x 10~ 3 rad/sec. 

• The reference density is taken to be po = 10 3 kg/m 3 . 

We then calculate the nonlinearity parameter ([3]) and the nonlinear time scale ([2]). To do 
so we substitute the Garrett and Mu nk spectrum (|2"1[) into the kinetic equati on with broaden- 
i ng (fl8l) . For matrix elements we use lMuller and Olbersl (119751 ). Eq. (jSJ), and lLvov and Tabak 
( 12004 ). Eq. ffT2|) . We also use the dispersion relation of internal waves, (j3J) for the isopycnal 
Hamiltonian, and ([6]) for Lagrangian coordinates. We use two values of 5 in (T221 : 5 = 10~ 2 
and 5 = 10~ 3 . We therefore make four calculations: 



Run I iLvov and Tabakl fl2004h with 5 = 10" 3 
Run II iMiiller and Olbersl fll975h with 5 = 10~ 3 



Run HI ILvov and Tabakl (J200J) with 5 = 10~ 2 
Run IV lMiiller and Olbersl dl975h with 5 = 10~ 2 



Results appear in Figs. [5] and El 

For Run 1 the nonlinearity parameter is uniformly small, smaller than 10 _1 . Such value 
of the nonlinearity parameter indicates that the kinetic equation is a self-consistent approach 
for the Garrett and Munk Spectrum. Increasing values of the nonlinearity parameter are 
noted with increasing vertical wavenumbers. This is consistent with intuition that we have 
about such systems. The nonlinear time scale is of the order of one hundred wave periods 
at low vertical wavenumber and of order ten wave periods at high vertical wavenumber. We 
also define a "zero curve" - It is the locus of wavenumber-frequency where the nonlinearity 
parameter and time-derivative of waveaction is exactly zero. The zero curve clearly delineates 
a pattern of energy gain for frequencies / < u> < 2f, energy loss for frequencies 2f < uj < 5f 
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and energy gain for frequencies 5/ < lu < N. This seems to be a characteristic pattern that 
appears in our calculations. Note that the zero curves are nearly independent of vertical 
wavenu mber. 

The iMiiller and Olbersl (119751 ). matrix element ([8]), Run II (5 = 10~ 3 ) results are quali- 
tatively similar to Run I. Factor of 2-3 faster decorrelation times and levels of nonlinear ity 
are noted in the high-frequency and high-wavenumber part of the spectrum. 

Therefore we conclude that when near-resonant interactions are included, the transfer 
rates are representation dependent. Furthermore, Lagrangian approaches predict higher 
level of nonlinearity. 

To investigate in more details results of near-resonant interactions, we perform numerical 
calculations for 5 = 10 -2 . Results for the canonical Hamiltonian formulation in isopycnal 
coordinates Run III are nearly identical to those with 5 = 10~ 3 . Results for the Lagrangian 
coordinate representation are both quantitatively and qualitatively different. The Lagrangian 
coordinate formulation ([8]) now predicts 0(1) nonlinearity for high frequencies, while the 
isopycnal coordinate formulation still returns nonlinearity parameter and much slower decor- 
relation time estimates. The zero curves for the Lagrangian coordinate representation are 
no longer simple functions of frequency at this higher level of nonlinearity. The zero curves 
in the isopycnal coordinate system are relatively independent of 5. 

To investigate the differences between approaches and the sensitivity of our results to the 
value of 5, in more details, we plot in Fig. [7] the differences of the nonlinearity parameter 
for these runs. In particular, we calculate the differences between Run I and Run II, Run I 
and Run III, and finally between Run II and Run IV. Differences associated with increased 
resonance broadening are minimal, 10~ 3 or smaller, for the isopycnal Hamiltonian. As the 
nonlinearity parameter estimates are representation dependent, differences between isopycnal 
coordinate and Lagrangian coordinate representations are much larger and increase with 
increasing 5. 

We have found that transports for the canonical Hamiltonian representation are not too 
sensitive to near-resonant interactions. We have also found in Section H] that all approaches 
are equivalent on the resonant manifold. We therefore conclude that all approaches will 
converge to Hamiltonian one as delta decreases. We have not undertaken such calculations 
as such small values of 5 would require significant modifications to our numerical algorithm. 
Note that the Fig. El especially Runs I and II, bear a strong resemblance to Fig. 4 of 



Olbersl (119761 ). These figures contain two positive and one negative lobe with similar bound- 



aries separati ng these regions, consistent with the characteristic pattern mentioned above. 



Olbersl (119761 ) does not make the hydrostatic approximation, used the GM75 model as the 



basis of his evaluations and is constrained to the resonance manifold. We have made the 
hydrostatic approximation, base our evaluations on the GM76 m odel and have included reso- 
nanc e broadening. Similarities are also apparent with Fig. 12 of (jMcComas and Brethertonl . 



19771 ) and Fig.s 10 and 11 of McComas and Miiller (1981). In those resonant evaluations 
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using the GM76 model, the hydrostatic approximation was invoked and interactions with 
frequencies greater than iV/3 were excluded. The major difference is that the zero line 
separating the positive and negative lobe at high frequencies has moved to 10/. 



6 Discussion 



In this paper we have review different approaches for wave-wave interactions that have 
been presented in the literature in the last t hree decades . Nam e ly, we have concent r ate o n 
the approaches of IMuller and Olbersl dl975l ): IVoronovichl (jl979h : ICaillol and Zeitlinl (j2f)00h : 
Lvov and Tabakl (120011 . 120041 ). In the absence of background rotation, we demonstrate that 
these four approaches produce equivalent results on the resonant manifold. 

This statement is not trivial given the different assumptions and coordinate systems 
that have been used for the derivation of the various kinetic equations. It points to an 
internal consistency on the resonant manifold that we still do not completely understand 
and appreciate. 



This result is l ess su rprising for the canonical Hamiltonian approaches (IVoronovichl . 11979 



Lvov and Tabakl . 120011 ). A canonical Hamiltonian representation is the gold standard of 
wave turbulence. It guarantees that the symmetries and hence conservation principles of the 
original equ ation set in the spatia l/temporal domains have been preserved in the spectral 
domain, (e.g. Zakharov et all Il992l ). Thus, if Voronovich's Clebsh variable representation in 
Eulerian coordinates and the Lvov and Tabak isopycnal Hamiltonian describe the same phys- 
ical system, then there is a canonical transformation that connects these two Hamiltonians. 
It is well known that such a canonical transformation reduces to the identity transformation 
on the resonant manifold. To prove this statement one constructs a near-identical canoni- 
cal transformat i on, wh ich is applicable for weakly nonlinear systems (See Appendix A3 in 
Zakharov et al. (fl99~l n. The Hamiltonian on the resonant manifold is invariant under a 
canonical near-identity transformation. 

That is why Voronovich's matrix elements (fTOl) look identical to the interaction matrix 
element in isopycnal coordinates ([TBI on the resonant manifold. 



We can argue that, w hile the other two matrix elements (ICaillol and Zeitlinl (120001 ) and 
Miiller and Olbersl (119751 )) are not in a canonical Hamiltonian formulation, they nevertheless 
do describe the same physical system. Consequently, they also can be approximated by a 
certain Hamiltonian structure, at least for small nonl i nearit ies. This is explicitly the case for 
the no ncanonical Hamiltonian o f (IMuller and Olbersl . Il975l ). It appears to be implicitly true 
of the ICaillol and Zeitlinl d2000h non-Hamiltonian kinetic equation. Therefore equivalence 
of the scattering matrix element on the resonant manifold is an intuitive, yet not trivial 
statement. 

On the other hand, it is also intuitive that there will be coordinate representation de- 
pendent differences. It is a robust observational fact that Eulerian frequency spectra at 
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high vertical wavenumber are contaminated by vertical Doppler shifting: near-inertial fre- 
quency energy is Doppler shifted to higher frequency at approximately the same verti- 
cal wavelength. Use of an isopycnal coordinate system considerably reduces this artifact 
(ISherman and Pinkell. Il99lh. Further differences are anticipated in a fully Lagrangian co- 
ordinate system (jPinkel . 2008 ). Thus differences in the approaches may represent physical 
effects and what is a stationary state in one coordinate system may not be a stationary 
state in another. In particular, differences may represent the effects of resonance broadening 
associated with Doppler shifting. 

We also demonstrate that the isopycnal and Lagrangian coordinate system approaches 
predict qualitatively different results with the inclusion of the near-resonant interactions and 
background rotation. The canonical Hamiltonian isopycnal formalism is insensitive to off- 
resonant interactions: Broadening the resonance width by an order of magnitude does not 
create significant differences in the nonlinearity parameter. The noncanonical Lagrangian 
coordinate representation is, in contrast, quite sensitive to these changes. 

As explained above, the Hamiltonian on the resonant manifold is invariant under near- 
identity canonical transformations. The kinetic equation describes the spectral transfers 
associated with the cubic terms of the Hamiltonian and conserves the energy associated 
with quadratic terms of the Hamiltonian. The kinetic equation should therefore return 
representation independent results on the resonant manifold. This statement is no longer 
true for near-resonant interactions. 

Indeed, since the structure of the Hamiltonian may be altered off the resonant manifold 
by a near-identity canonical transformations, one should anticipate representation dependent 
differences in spectral energy transfer when near-resonant interactions are included. Such 
differences become more and more significant as nonlinearity increases and cubic parts of 
the Hamiltonian become increasingly large. 

We would like to s uggest that the differences between IMiiller and Olbersl (119751 ) and 



Lvov and Tabaki (120041 ) off the resonant manifold represent physical effects. However, an 



issue with extant Lagrangian coordinate representations is that they require a small ampli- 
tude assumption that represents an unconstrained approximation whose domain of va lidity 



vis-a-vis the weak interaction approximation is not well defined, (iMiiller et all Il986l ). On 



the basis of estimates of how horizontal Doppler shifting contributes to isopycnal spectra, 
we would anticipate that the Lagrangian coordinate stationary state would have typically 
steeper spectral slopes in the frequency domain than frequency spectra in isopycnal coor- 
dinates. The results presented here indicate that resonance broadening will quickly whiten 
the high frequency Lagrangian coordinate spectrum, in direct contradiction to our intuition 
regarding physical effects. 

In this paper we have shown that while on the resonant manifold (i.e. for weakly nonlinear 
interactions) all approaches we considered do agree, inclusion of the near-resonant interac- 
tions (for stronger nonlinearities) should be done with care. Results with near resonant 
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interactions are representation dependent. This observations warrants further study. 
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Figure 1: Matrix elements \V£ 



PUP2 | 



15bt ! 2 given by the solution (pb]). upper left: l ^^jj^l 



according to iMiiller and Olberd (119751 ). upper rig ht: l^,p 2 |Y5|^| 2 ac cording to IVoronovich 
f ll979f ). bottom left: l^p 1 ;P2 fe^| 2 accordi ng to ICaillol and Zeitlinl (l2000h . bottom right: 
l^pi.P2^5b|| 2 according to lLvov and Tabakl (1200 ll ). 
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Figure 2: Matrix elements iV^p JigTj]! 2 given by the solution (I16al) . upper left: | V^ p rj-g^|| 2 
according to iMiiller and Pipers! (I1975I ). upper ri ght: Ij^pjjg^l 2 a ccording to IVoronovich 
fll979f ). bottom left: I ^Jfal 2 accord ing to ICaillol and Zeitlinl fl2000h . bottom right: 
\ V plpMa[\ 2 accordin g to lLvov and Tabakl fl200ll ). 
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Figure 3: Matrix elements | ^^ rjg^ | 2 given by the solution ( l!6bl) . upper left: l ^^pfig^ll 2 
according to iMiiller and Olberd (I1975I ). upper ri ght: l^^plxgbjj 2 a ccording to IVoronovich 
fll979f ). bottom left: I V ^pf^ 1 2 accord ing to ICaillol and Zeitlinl fl2000f ). bottom right: 
iV^pfebll 2 according to lLvov and Tabakl (l200lh . 
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Figure 4: upper: Matrix elements l^ P2ES | 2 given by the solution (j!5ap . middle: Matrix 
elements |V^ P2ID | 2 given by the solution (Il5bl) . bottom: Matrix elements l^^ppgj 2 given 
by the solution (116aj) . which gives PSI as |fejj ^ (e — > 0). The matrix elements here are 
shown as functions of e such that (|fei|, |fc 2 |) ^(e, e/3 + l)|fe|. All four versions of the Matrix 
elements are plotted here: the appearance of a single line in each figure panel testifies to the 
similarity of the elements on the resonant manifold. 
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Figure 5: Nonlinear ity parameter ([3]) for the Garrett and Munk spectrum (f24j) calculated 
via f[T8|). The upper fig ures represent the value of nonlinearity parameter calculated using 
Lvov and Tabakl (120041 1 . equation ( !T2l with 5 = 10~ 3 , Run I (upper left) and 5 = 10 -2 Run 



III (upper rig ht). The bottom two pi ctures represent the value of nonlinearity parameter 
calculated via lMuller and Olberd fll975h . Q3Q with 5 = 10" 3 Run II (bottom left) and S = lO" 2 
Run IV (bottom right). 
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Figure 6: Nonlinear time ([2]) for the Garrett and Munk spectrum (12411 calculated via 
( fT8l) . The uppe r figur es represent the value of nonlinearity parameter calculated using 
Lvov and Tabakl (l2004f ). equation QE) with 5 = 10~ 3 , Run I (upper left) and 5 = 1(T 2 
Run III (upper rig ht). The bottom two pictu res represent the value of nonlinearity param- 
eter calculated via Muller and Oibersl (119751 ). flS} with 5 = KT 3 Run II (bottom left) and 
5 = 10~ 2 Run IV (bottom right). On this bottom right figure white region to the left of 
the 0.1 contour corresponds to extremely fast time scales, faster then 0.1 of a day. On these 
figures, uj in cpd, m in cycle/m, and nonlinear time r NL is measured in days. 
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Figure 7: Differences betwee n nonl inearity parameter 

3 



calculated via Miiller and Olbersl 
(ll975l ) and iLvov and Tabakl (l2004f) wi th 5 = 10~ 3 , i.e. between Run I and Run II, (a), 
calculated with lLvov and Tabakl (I2004T) with S = 10~ 2 and S = 10 -3 , i.e. the difference Run 
I and Run III (b), and finally between iMiiller and Olbersl (Il975l ) with 5 = 10 -2 and 5 = 10 -3 , 
i.e. between Run II and Run IV (c). 
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